Quantum site percolation on triangular lattice and the integer quantum Hall effect 
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Generic classical electron motion in a strong perpendicular magnetic field and random potential 
reduces to the bond percolation on a square lattice. Here we point out that for certain smooth 
2D potentials with 120° rotational symmetry this problem reduces to the site percolation on a 
triangular lattice. We use this observation to develop an approximate analytical description of 
the integer quantum Hall transition. For this purpose we devise a quantum generalization of the 
real-space renormalization group (RG) treatment of the site percolation on the triangular lattice. 
In quantum case, the RG transformation describes the evolution of the distribution of the 3x3 
scattering matrices at the sites. We find the fixed point of this distribution and use it to determine 
the critical exponent, v, for which we find the value v « 2.3 2.76. The RG step involves only a 
single Hikami box, and thus can serve as a minimal RG description of the quantum Hall transition. 

PACS numbers: 72.15.Rn, 73.20.Fz, 73.43.-f 
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I. INTRODUCTION 

Network-model formulation of the Anderson localiza- 
tion problem was first introduced in Ref. d The key 
observation made in Ref. IT] is that a complex motion 
of electron in disorder potential can be reduced to the 
motion along the links of the network (in both direc- 
tions) with disorder incorporated via random phases of 
scattering from the nodes of the network. Since then, 
the network-model approach became a powerful tool for 
numerical studies of disordered systems. In these stud- 
ies the randomness is incorporated into the phases ac- 
cumulated along the network links. A great advantage 
of the network-model approach to localization is that 
it can be conveniently applied to various universality 
classes^ of disorder. In order to capture the specifics 
of a given class, one has to impose an appropriate sym- 
metry requirements on the random phases on the links, 
which fixes the form of 5'-matrix at the nodes. This is 
achieved by introduction of an internal space associated 
with the link of the network model, which possesses a 
desired symmetry. For example,"^ with two-component 
links (one component for one spin projection) for each 
direction of propagation the requirements of unitarity 
and time-reversal symmetry allowed to reveal a delocal- 
ization transition expected in two-dimensional systems 
with spin-orbit scattering. Comprehensive review of the 
results obtained with the help of the network model is 
given in Ref. [5|. In addition to establishing the exis- 
tence of localization transitions in different classes, nu- 
merical simulations of the network models with the help 
of transfer-matrix method yields quantitative character- 
istics of the critical point. These characteristics include 
critical exponent, critical level statistics, and critical con- 
ductance distribution. Corresponding references can be 
found in the review Ref. [Sj. Note that, with the exception 
of Ref. i where a particular two-channel network model 
was considered for arbitrary graph, the underlying net- 
work in all previous studies was a square lattice. 

Especially convenient for modeling with a network is 
the chiral motion of a 2D electron in a disorder poten- 



tial and a strong perpendicular magnetic field. This is 
because the corresponding network is directed. Directed 
character of the network with chiral scattering at the 
nodes allowed Chalker and Coddington^ to demonstrate 
conclusively that there is only a single delocalized state 
per Landau level. This finding is of a great importance, 
since it is the underlying reason for sharp conductivity 
peaks at the quantum Hall transition. 

On physical grounds, the seminal Chalker-Coddington 
modelS of the integer quantum Hall transition can be 
introduced in a natural way with the help of a two- 
dimensional potential 



U {x, y) — cos(7ra;) cos(7ry). 



(1) 



In this potential, the equipotcntial lines U (x, y) = Q 
form a square lattice. For any nonzero e equipotentials 
U{x, y) = s are closed. For positive e these equipoten- 
tials encircle the maxima (x, y) = (2to, 2n) [and also 
(2to + l,2n + 1)] of U{x,y), while for e < equipo- 
tentials encircle the minima {x,y) = (2m, 2n -I- 1) [and 
(2m + 1, 2n)] of U{x, y). In a strong perpendicular mag- 
netic field 2D electron drifts along equipotentials. Then 
reconfiguration of equipotentials at e = 0, as illustrated 
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FIG. 1: (Color online) (a) Evolution of equipotentials near 
the saddle point of potential Eq. ^ as the energy e crosses 
over from e < to £ > 0; (b) evolution of equipotentials near 
the nodes of potential Eq. with e. Three equipotential 
lines touch at e = 0. 
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FIG. 2: (Color online) Equipotontial lines e = of the poten- 
tial Eq. Q constitute a triangular lattice. Minima of W{x, y) 
are in the centers of shaded triangles. Blue lines constitute a 
triangular lattice of supersites with lattice constant, \/3. Red 
lines: hexagonal lattice with sites in the centers of triangles. 
Each site of a triangular lattice is surrounded by a hexagon. 



in Fig. [T^, manifests the change in the character of mo- 
tion. 

Chalker and Coddington^ captured the quantum char- 
acter of motion in U (x, y) by assigning to the saddle 



points at {x,n,yn) = [m 
tial behaves as 



U{x - Xm,y - yn) ~ 

(-l)'" + "7r^(x - Xm){y - yn) 

a scattering matrix 
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Here p and if are the polar coordinates with origin at 
(xm ,yn), and e is the dimensionless energy. For a realis- 
tic potential with magnitude, J/q, and correlation length, 
d, dimensionless energy, e, is the physical energy mea- 
sured in the units of the width, F ~ UqP /(P, where, 
is the magnetic length. Chalker and Coddington demon- 
strated that, in order to account for a smooth disorder, 
it is sufficient to assume that the phases, acquired by a 
drifting electron between the saddle points, are random. 

In this paper we show that the description of the quan- 
tum Hall transition can also be obtained based on the 
potential 



W{x, y) = ^ Vij — nei — me2). 



(4) 



which has a triangular symmetry, i.e., ei = (1,0) and 
62 = (1/2, \/3/2) are the basis unit vectors of a triangular 
lattice; the function V is defined inside a "black half" of 
the rhombus-shape unit cell, see Fig. ^ in the following 
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Here p and ip are the polar coordinates with respect to the 
origin at (ei +62)/^. The form of potential in the "white 
half" of a unit cell is given by Eq. ^ upon replacement 
V{p, ip) —V{p, p — 7r/3) and shifting the origin to the 
center of the white triangle. 

Equipotentials of Eq. (O evolve, as e passes through 
0, in a fashion qualitatively different from the case of 
quadratic symmetry. As shown in Fig. [T^, in the case 
of quadratic symmetry, black regions (minima) get con- 
nected, while adjacent white regions (maxima) get dis- 
connected. By contrast, the behavior of the potential 
W{x,y) near the nodes at r,„_„ = (m — l/3)ei + {n ~ 
l/3)e2 is given by 



Wip, p) = -p" J2 sin [3(2g + l)p] 



g=0 



i) where the poten- where 



2\/3 



[(2g + l)(9(2g- 
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(6) 



Corre- 



sponding evolution of equipotentials is illustrated in 
Fig- [lb- We see that, as e crosses zero, three black regions 
get joined at e = simultaneously. This suggests that 
quantum mechanical description of motion in the poten- 
tial W{x^y) requires, in addition to random phases on 
the links, introducing a 3 x 3 scattering matrix at each 
node. Below we argue that the form of this matrix is 



1(1 
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(i-e)\ 



) 1(1 + ^) 



(7) 



(1-e) |(l+e)/ 



With matrices Eq. ([7]) in the nodes, the correspond- 
ing network model is shown in Fig. [3] In this net- 
work the phases on the links are random, as in Chalker- 
Coddington model, while all S'-matrices in the nodes, 
Eq. ([7]), are the same. As we demonstrate in present 
paper, this model can be treated numerically using the 
same MacKinnon-Kramer finite-size scaling algorithm^ 
that was employed in Ref. H (for subsequent numerical 
studies of the Chalker-Coddington model see Refs. 
ll0!lll|il2lMJ4,1516^17,18^19,20 and review articles^). 

Chalker-Coddington model can be viewed as quantum 
version of the classical bond percolation. Establishing 
one to one correspondence between the classical bonds 
and the links is possible due to directed character of the 
chiral CC network. On the other hand, it was demon- 
strated in Ref. [13 that a simple real-space renormaliza- 
tion group (RG) procedure, based on decimation, leads 
to the closed equation for the classical bond percolation 
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FIG. 3: Network model on a triangular lattice is illustrated. 
Three links enter a node and three links exit each node; the 
nodes with S-matrices Eq. (O are depicted as triangles. 



threshold. This procedure reproduces the exact thresh- 
old value and yields a very accurate estimate for the crit- 
ical exponent. In Refs. 23,24,25,26,23 the classical RG 
procedure^^ was generalized to the quantum bond per- 
colation. It was shown that corresponding integral RG 
equation yields, in addition to the accurate value of the 
quantum critical exponent, a very accurate distribution 
of the critical conductance. A real-space RG procedure 
for a site percolation was proposed in the same paper Ref. 
[2^ . This procedure is simpler than for the bond perco- 
lation and yields more accurate results. In the present 
paper this procedure is extended to the quantum case, 
where it describes the critical behavior of the directed 
network Fig. [31 As in the classical case, the quantum RG 
analysis of on the triangular lattice is much easier than 
the quantum RG analysis of the Chalker-Coddington on 
the square lattice. This analysis is presented in Sections 
II-V. 



II. RG PROCEDURE FOR CLASSICAL 
PERCOLATION ON A TRIANGULAR 
NETWORK 

In order to develop the RG description, it is necessary 
to incorporate disorder into the ^'-matrices. The reason 
is that at each RG-step three ^-matrices, connected by 
the links, are combined into one swper-S'-matrix. Then 
the randomness in phases translates into the randomness 
of the elements of super- S-matrbi. 

A natural way to devise the quantum RG description 
is to start from a classical problem of electron drift along 
the sides of triangles. Fig. O 



A. Classical drift on triangular lattice 

Assume that the potential W{x, y) is perturbed near 
the nodes as 

W{p, = Vb - cop2 sin(3^) - cip2 sin(9(p) , (8) 



where the random shift, Vq, is much smaller than 1, but 
much bigger than the "quantum" energy width, F. De- 
pending on the shift, an electron drifting along equipo- 
tential W ^ e towards a node, turns either to the left 
(Vo < e) or to the right (Vq > e), see Fig. [TJ). We can 
conventionally call the nodes with Vq < e and Vq > e 
as "black" and "white" lattice sites, respectively. It is 
obvious that when the average, (Vq), is zero, the perco- 
lation threshold in the potential Eq. ([8]) is e = 0. In the 
language of sites, the same is to say that 50% of sites are 
black at e = 0. It is, in fact, well-known^§. that the ex- 
act threshold for site percolation on a triangular lattice is 
50%. However, it is much less obvious that there is a com- 
plete equivalence between electron drift in potential Eq. 
([5]) and the site percolation on a triangular lattice. Su- 
perficially, this can be seen from the fact that when two 
neighboring sites are black, they are connected via the 
black region. Still, a rigorous proof requires additional 
steps, in particular, introducing auxiliary hexagons, see 
Fig. O This proof is given in the Appendix A. 



B. Classical RG scheme 

The above mapping on the percolation problem allows 
one to employ the RG approach to the site percolation on 
triangular lattice put forward in Ref. |2^ . This procedure 
is much simpler than the RG for the bond percolation on 
the square lattice, proposed in the same paper. Note in 
passing, that the bond percolation on the square lattice 
is the classical limit of the Chalker-Coddington modcl^'^.. 

As shown in Fig. [21 at each RG step2^ the lattice 
constant increases by a factor of \/3. A site of a rescaled 
lattice, a supersite, is either black or white depending on 
the colors of the three constituting sites: if either all three 
or only two out of three constituting sites are black, then 
the supersite is black. Otherwise, the supersite is white. 
Quantitatively, the probability p' for a supersite to be 
black is expressed via the corresponding probability for 
the original site as 



P 



R{p)^p^ + 3p^il~p). 



(9) 



Fixed point, p' = p = Pc, of the transformation Eq. ^ 
reproduces the exact result Pc = 1/2. Critical exponent 
is determined by the condition that the correlation ra- 
dius, ^ = {p ~ Pc)~'^ , on the original lattice is equal to 
the correlation radius ^/S {p' — pc)^'^ on the renormalized 
lattice, i.e. 



In 



In {V3) _ In (%/3) 
^ In (3/2) ■ 



\ ^ / p—Pc 



(10) 



Eq. PH)) yields ly ~ 1.354, which differs from the exact 
value, 1^ — 4/3, by only 1.6%. 

The rationale behind the transformation Eq. ([9]) is 
that the supersite is located in the center of the black 
cell in Fig. [21 Then the color of the supersite reflects 
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the "percolation ability" of this black triangular cell, so 
that, even if one of the nodes, constituting the vertices of 
the triangular cell is white, the cell still percolates "over 
black". 



C. RG in the language of potential shifts 

At this point we make an observation that the above 
RG procedure can be reformulated in the language of 
potential W{x,y) with random shifts in the nodes, Vq. 
Namely, for Vg^, Vq, and V^^ being the shifts at the nodes 
constituting a supernode, the shift Vq of the supernode 
is defined as 



V^ = Mid{Vo\Vo'X}, 



(11) 



where Mid stands for Vq which is smaller than maximal 
but larger the minimal out of the three numbers. With Vq 
defined by Eq. (fTT|) . the RG equation Eq. Q describes 
the evolution of probability that the shift exceeds e. 

The importance of the above observation is that re- 
formulation of classical RG procedure in terms of poten- 
tial shifts opens a possibility to capture the quantum- 
mechanical motion in the random potential. A prescrip- 
tion how to extend classical description to the quantum 
case is^'^: the algorithm Eq. ITTI) should be cast in the 
form of a scattering problem. 



III. REFORMULATION IN TERMS OF 
CLASSICAL TRANSMISSION 

We identify the scattering object as a point where three 
equipotentials come close as shown in Fig. 21 Incident 
electron, ii, either proceeds along the same equipoten- 
tial into oi (reflection) or switches equipotentials and 
proceeds along 02 . Retention of equipotential (reflection) 
corresponds to positive Vq — £ in the vertex, encircled in 
Fig. m In terms of scattering matrix 











03 / 





(12) 



the same simple notion can be reformulated as follows. 
For positive Vq — e, the matrix S = 5*+ is a unit matrix, 
while for negative Vq — e the matrix S ^ S- has the form 



Superscattering object consists of three scattering ob- 
jects, and is shown in the same figure. 

Now we reformulate Eqs. ([9]) and pT|) in yet another 
language of ^'-matrices. Namely, the ^-matrix of a su- 
perscattering object is expressed via S'-matrices of con- 
stituting scattering objects upon reducing the number 




of legs from 12 to 6. We emphasize that this reduction 
can be carried out in two distinct ways, as illustrated 
in Fig. O The first way is to perform contractions as 
02 <-> ^2, 04 <-> «4, and og ^ ig. The second variant of 
contractions is 02 <-*■ ii, 04 13, and og 15. Now 
it is straightforward to check that RG transformation 
Eq. (jlip corresponds to the following rule for 5-matrix 
of the superscattering object, S. If the 5-matrices of ei- 
ther all three or only two of constituting objects are S+, 
then S — 5*+. In all other realizations, when at least 
two of constituting objects have the matrix we have 
S = 

It is important to note that the above rule for S applies 
independently of the way in which the contractions in 
Fig. O are performed. This is not the case in the quantum 
version to which we now turn. 



IV. QUANTUM GENERALIZATION 

Quantum 5'-matrix of the scattering object differs from 
5+ and S- in two respects. Firstly, at the points of close 
contact between each pair of equipotentials electron can 
switch equipotential even for positive Vq — e, when it is 
forbidden classically. Corresponding classically forbidden 
transitions are illustrated in Fig. 2] with red dashes. Sec- 
ondly, upon travelling between two subsequent points of 
close contact, electron accumulates the Aharonov-Bohm 
phase, (pi. For example, the phase 0i is accumulated in 
course of drift between ii and oi. These phases are irrele- 
vant in the classical limit when the reflection amplitudes, 
yi, are either or 1. However, for intermediate < < 1 
the amplitude for an electron to execute a close contour 
around the center in Fig. [4] (following the red dashes in 
the clockwise direction) is finite. As a result, 4>i enter 
into quantum scattering matrix. Explicit form of S in 
terms of < j/i < 1 and can be obtained by solv- 
ing three pairs of linear equations, describing quantum 
scattering at each of three points of the close contact of 
equipotentials. We have 



S = 



-1 



1 — xiX2X3e^'^ 



(13) 



/ yiy3e"*i a;i?/22/3e'<*i+*=' xiX2e'^ - X3 \ 

X2X3e''* - xi yme"^^ yiX2y3e'<*=^+*^' 
V yiy2a;3e'<*i+*3' xiX3e^-x2 ^22/36'*=' / 



where Xj = 
tudes, and 



1 — stand for the transmission ampli- 



Tp ^ (pi + (f)2 + (f>3 



(14) 



is the net phase accumulated along the closed contour. 
It is easy to check that the matrix Eq. (|13p is unitary. 
It is also straightforward to verify that in the classical 
limits, when all Xj = or = 1, Eq. (|13p correctly 





FIG. 4: (Color online) Left: superscattering object is shown 
with a big red dot, while constituting scattering objects are 
shown with small red dots. Right: detailed structure of the 
scattering object; black dots show the points of the close con- 
tact of equipotentials. Electron incident along equipotential, 
i, can either proceed along the same equipotential, a, or swich 



equipotential with probability 
illustrated with red dashes. 



a; I . The latter processes are 





a 



FIG. 5: (Color online) Two contributions to the kernel Eq. 
(|19p originating from two variants of contraction of equipo- 
tentials are illustrated; (a) corresponds to i", Eq. (|16|) . while 
(b) corresponds to x'', Eq. (flT)) . 



B. Quantum RG equation 



reproduces S-i- and respectively. Detailed derivation 
of the form Eq. of the scattering matrix is presented 
in Appendix B. 



The form of the matrix Sa 



From Eq. ([13)) we can establish the form of the scatter- 
ing matrix St,{£), Eq. (O, with the help of the following 
duality argument. Due to triangular symmetry of the 
potential sin(3<y9) we have xi = X2 = X3 = x. Consider 
now the transition point, e 0. At this point the prob- 
ability for electron incident along, say, ii (see Fig. S]) to 
be deflected to the left (along oi) is equal to the proba- 
bility to be deflected to the right (along 02). Less trivial 
is that the phase, ip, must be zero at e = 0. This is the 
consequence of the fact that the scattering problems for 
electron with energy e and — e are equivalent if we change 
the drift direction from clockwise to anti-clockwise; this 
change implies also the change in the sign of ip. 

With V = 0, the condition, (1 — x^)^ = {x^ — x)^, of 
"equal deflection" to the left and to the right, yields a 
single physical root x — —1/2. Substituting it back in 
the matrix Eq. (|13p reduces it to the scattering matrix 
Eq. ([7]) with £ = 0. Then for probabilities ii oi 
and ii 02 we get 4/9, while the probability ii —^ 03 
is 1/9. Including small finite e can be also performed 
with the help of the duality argument, namely that the 
probability of deflection to the left with energy e is equal 
to the probability of deflection to the right with energy 
— e. On the other hand, the change of probability ii — > 03 
with e is (X e^. 

The role of the matrix 5a is central to the transfer- 
matrix treatment of the network model; see Appendix 
C. 



We now turn to the quantum RG procedure. In con- 
trast to evolution of probability, p, upon rescaling of the 
lattice constant in the classical case, this procedure is 
formulated in terms of evolution of the distribution func- 
tion, Q{x), of the absolute values of the transmission am- 
plitudes, Xi. Thus the quantum generalization of Eq. ^ 
at the step n is the following recurrence relation 



Q7i+l{x) = T {Qn{x)} 



(15) 



J ^^Y\dXjQn{Xj)^ K{x,Xi,X2,X3). 

The kernel, K, represents the conditional probability 
that, after performing the contractions, the transmission 
coefficient v ~^ oi in Fig. 3] is equal to x, provided that 
the constituting transmission coefficients are xi, X2, and 
illustrated in Fig. [D In analytical evalu- 
ation of the dependence x{xi, X2, X3) it is important to 
take into account that this dependence is different for two 
variants of contractions. For a variant 02 i2, shown in 
Fig. [5^, the transmission coefficient is given by 



(16) 



1 - 



(1 - xiX2X3e^-^)e'v + ^(^i^xl){l-xl) 



where ip is the phase along the contour (02,^2) , which 
is now closed. Correspondingly, for the second kind of 
closing 02 ii. Fig. Os, the transmission coefficient has 
the form 



(17) 



1 - 



i^,-e^^)Vil-xi){l-xi) 



(1 — xiX2X^e^'^)e^^ + X2X^e 



lip 



Xl 



where f stands for the phase along (02,^1). For details 
of derivation of Eqs. ^17} see Appendix B. 
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FIG. 6: (Color online) Left: Kernel K{x,xi,X2,X3) is plotted 
from Eq. (|19|l for three sets {xi,X2,X3). Right: approximate 
analytical solution (green line) of the RG equation Eq. (I21|l . 
given by Eq. (1^ . is plotted together with r.h.s. of Eq. (1^ . 
which is T {Q{x)} (black line). 



Relations Eqs. (fT7|) define also the dependencies 

x1'^{xi,X2,X3) and X2''''{xi,X2,X3). It is important that 
the central phase, ip, is common for all three dependen- 
cies {ij({xi},(pj, V')}j=i,2,3 and {^^({xi}, (p^ , i/')}j=i.2.3- 
A crucial step in extending classical RG procedure to 
the quantum case is to follow the " classical" prescription 
to choose a middle out of three coefficients 

i'^({xJ,{^,},V^)= (18) 

Mid(x?({xJ,V3l,V'), X^{{x^},ip2,1p), S3({X^},(^3,V')), 

and the same for a;''({a;i}, {iy9i}, ■(/;). Eq. (|18p is a quantum 
generalization of the classical Eq. pT|) . Note that selec- 
tion of middle value in Eq. (|18p is performed for given 
values of random phases, ipi, ip2, ifs, and "0- Within RG 
procedure different phases are uncorrelated, and in eval- 
uation of the kernel we average over each of four phases 
independently. Finally, taking into account that contrac- 
tions a and b are statistically equivalent, the expression 
for K{x, {xi}) acquires the form 

K{X, {X,}) = Usix - i'^({x.}, V))) 



6{x~S: {{xi},{ip^},ip)) 



(19) 



Quantum delocalization corresponds to the fixed point 
of the transformation Eq. ([TS]). It is found in the next 
Section. 



NUMERICAL RESULTS FOR FIXED POINT 
AND CRITICAL EXPONENT 




FIG. 7: (Color online) T,,, = T|q(2 - ^o)} is plotted for 
different zq. Inset shows the variance, x(''')i plotted from Eq. 



(|26p . The minima correspond to r 
zo — 1 and zq — —1, respectively. 



1.27 and t = 1.22 for 



for different sets xi, X2, X3. It is seen that the kernel sup- 
ports the attractive critical points Xi = X2 = X3 = and 
xi — X2 — X3 — I. Indeed, for values xi = X2 = x^ — 0.2, 
the kernel is centered at even smaller value x « 0.04, 
while for xi — X2 = X3 = 0.85 it is around bigger value 
X « 0.95. In both cases the kernel is narrow. This is 
because for "classical" transmission coefficients interfer- 
ence does not play a role, so that the phases drop out 
from Eqs. p^ . (|17p . However, for intermediate val- 
ues xi = 0.6, X2 = 0.65, X3 = 0.7 the kernel extends 
over entire interval < x < 1. It also exhibits peaks 
at a; = 0.03,0.61,0.65,0.8 and 0.83. The origin of these 
peaks is the anomalous contributions of phases (/Sq = 0, tt 
and ipQ — 0,TT to the kernel. Indeed, for these values of 
phases we have dx°"''/d(p = dx°''''/dip = 0, so that 



X'-'i^P^lp) = X°'''((/3o, Vo) + 



(V-Vo)' + 5Li"'' 



iip-VoY (20) 



It is easy to check that, when the quadratic form of 
second derivatives is negatively defined, the correspond- 
ing contribution to the kernel is oc In |a; — x°''''{ipQ,4'o)\- 
Therefore, the peaks in the kernel reflect the fact that 
two loops corresponding to phases (f) and ip are insuffi- 
cient for complete averaging. Still, the phase volume of 
the singular contributions is small, so that the fixed point 
of the transformation Eq. is not sensitive to these 

singularities. 



The kernel 



B. Fixed point 



The examples of the kernel, K(x), calculated using 
Mathematica from Eqs. (fTB|) - P^ . are plotted in Fig. [H] 



The fixed point, Q{x), of the quantum RG equation 
Eq. (jl5p . satisfies the following nonlinear integral equa- 
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tion 



Qix) = T{Qix)} 



Y]_dxjQ{xj) \K{x,xi,X2,X3), 



(21) 



which we have solved using Mathematica. Starting from 
initial distribution, Q{x) = 1, and performing analytical 
fit at each step of successive approximations, the follow- 
ing expression for the fixed point was obtained after the 
fourth step 



0.9 



0.9 



1 + 103x2 1 + 5- 102(1 -a;)2 
+0.63x^ + 0.55(1 - xf + 0.6. 



(22) 



This expression is plotted in Fig. [5Jd with a green line. 
One can judge on the accuracy of the approximate solu- 
tion Eq. ((22)) by substituting it into the right-hand side 
of Eq. (PT|) . The result, black line in Fig. [5)3, is indeed 
very close to Eq. (j^ . 

The fixed point solution rises upon approaching x = 
and X = 1. This behavior is inherited from the classical 
percolation. Note that similar behavior was found in 
Refs. [23II25I for the fixed point of the quantum bond 
percolation on the square lattice. Direct comparison with 
conductance distribution P{G) found in Refs. [23II25I can 
be performed using the relation P{G) = Q{VG)/2VG. 
This comparison indicates that, while numerically the 
fixed point distributions are close, Eq. ([22]) favors small 
values of x. Qualitatively, this reflects the fact that at 
critical energy, e = 0, electron incident along ii (see Fig. 
HI) is more likely to proceed along ai rather than switch to 
02- The asymmetry is seen more clearly if one interprets 
Q{x) in terms of distribution, Q of heights of the effective 
saddle point. This height is determined by the relation 



l/(e^ 



1), so that 



2(e^ + 1)3/2 



([e^ + (23) 



The distribution Q{z) is shown in Fig. [7| with a black 
line. It has an asymmetry towards large z. 



C. Critical exponent 

To estimate the critical exponent, v, which governs 
the divergence of localization radius, ^(e) cx l/e", as 
a function of energy, e, we used the reasoning from 
Refs. 2 3,, 2 5, Electron with finite energy, e ^ F, "sees" 
the shifted distribution of the saddle point heights Q{z — 
zq), where <C 1 is proportional to e. The key step 
of the reasonin g23i25 jg that, upon the RG transforma- 
tion, the electron travels on the lattice with the lattice 
constant -\/3 and sees the shifted distribution of heights, 

T|(5(2 — zq)^ — Q{z — tzq), where r is some constant 



independent of Zq. After subsequent n RG steps this 
distribution evolves into Q{z — t"zo). When the shift 
accumulates to reach ^ 1, the electron becomes localized 
within the size of a unit cell of renormalized lattice. Then 
from the relations 

C(^o)(xi;, e=(\/3)", zot"^1, (24) 



we find 



lnV3 
Inr 



(25) 



This definition of is a quantum generalization 
of Eq. (fTO|) . In Fig. [7] the result of calculation 

T|q(2 — 2o)| for four zq = —1, —0.5, 0.5, 1 is shown. 

We see that for these zq the shape of T|q| is only 

slightly affected by the shift. The curves are approxi- 
mately equidistant, so that an estimate of r can be ob- 
tained simply from the horizontal separation « 0.4 be- 
tween the neighboring curves. This yields r « 1.25 and, 
correspondingly, v « 2.46. For more accurate estimate 
we studied the variance. 



X(r, ^o) = ^ [ T{q(z, - zo)} - Q{z^ - Tzo)] ' , (26) 



as a function of r for different values of the "energy 
shift" Zq. The sum Eq. (pS)) was taken over discrete set 
z, = -10 -f- O.Oli for i = 1, 2, .., 220. In Fig. [7] we plot 
the variance for zq = I and zq = —1. Both curves have 
pronounced minima at r = 1.27 and 1.22, respectively. 
This translates into the values of v — 2.3 and v = 2.76. 
Although these values are in good agreement with known 
value of the accuracy of the above estimate is limited. 
The limitation is due to the fact that for zq = ±1 the 

heights of maxima of the curves T|(3(2: — zo)|, shown 

in Fig. [3 differ from the height of Q{z). This devia- 
tion would not be a problem for smaller zq. However 
for zo = ±0.5 the variance becomes small, and its de- 
pendence on T becomes weak. Apparently, the variance 
Eq. ([26]) is affected by the wiggles at the top of the curves 

T|(5(2 — 2o)| much stronger for zq = ±0.5 than for 
Zq = ±1, which makes the evaluation of x for small zq 
ineffective. 



VI. CONCLUSION 

It is interesting to point out that, while the classical 
limit of the Chalker-Coddington model based on poten- 
tial Eq. ^ reduces to the bond percolation, similar form 
of potential Eq. ([8]) leads to the site percolation. The 
reason is the symmetry of corresponding potentials. As 
seen from Fig. ^ the hexagons surrounding the nodes 
of potential Eq. ([5]) have common sides. On the other 
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hand, the squares, drawn around the nodes of potential 
Eq. (21) share the vertices. 

Both the bond percolation on a square lattice and the 
site percolation on triangular lattice have pc — 1/2, which 
is insured by self-duality. The RG descriptions^- of both 
cases, having fixed point, = 1/2, effectively preserve 
this self-duality. As a result, the RG values for classical 
exponent come out close to = 4/3 in both RG schemes. 
In this paper we demonstrated that quantum extension 
of classical RG to the triangular lattice also yields the 
critical exponent close to the known value u = 2.33. 

The fact that the simple RG scheme, considered in 
the present paper, describes the quantum Hall transition 
so accurately has a deep underlying reason. Delocalized 
state in the quantum Hall transition emerges as a result 
of competition of two trends: (i) quantum interference 
processes that survive in the presence of magnetic field 
(Hikami boxes, Ref. [29l ) tend to localize electron, while 
(ii) classical Lorentz force, by causing electron drift, pre- 
vents it from repeating closed diffusive trajectories. Both 
trends are incorporated into our RG scheme. Obviously, 
chiral motion is the consequence of the Lorentz force. 
Hikami boxes, on the other hand, are represented in the 
RG step in the form of figure-eight trajectories, as illus- 
trated in Fig. m 

Note finally, that simplicity of the RG description 
proposed here suggest possibility to extend it to dif- 
ferent from qunatum Hall universality classes, see, e.g., 
Refs. [lillliilli 
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APPENDIX A 

Here we elaborate on the mapping of the problem of 
percolation over equipotential lines in the random poten- 
tial Eq. ([8]) and the conventional site percolation prob- 
lem. The easiest way to establish this mapping is to 
surround all sites of triangular lattice with hexagons, as 
shown in Fig. [51 If the site is occupied, then the hexagon 
is, say, black; if the site is vacant, it is white. The dis- 
tinctive property of the triangular lattice is that, when 
two hexagons touch, they automatically have a common 
side. Note in passing, that this is not the case for a square 
lattice, where two squares, drawn around the sites, may 
touch by sharing a vertex but not share a side. In the 
site-percolation problem, the bond between two neigh- 
boring sites conducts if both of them are occupied. The 
same is to say that conduction is possible between two 
touching hexagons, if they are both black. Percolation 
threshold corresponds to the portion of black hexagons 
when conduction over entire sample becomes possible. 



The fact that Pc = 1/2 is a direct consequence of geo- 
metrical arrangement of hexagons, due to which perco- 
lation over black hexagons rules out the percolation over 
white hexagons. Also, due to this arrangement, one of 
the colors always percolates. 

Consider now the potential Eq. Black sites are 

now those in which Vq < e. Configuration of equipo- 
tentials around this site is the rightmost of three shown 
in Fig. UId. Accordingly, configuration of equipotentials 
around the site with Vb > e is the leftmost of three shown 
in Fig. [T|d. Consider now two neighboring sites with 
Vq < e. Fig. [21 makes it apparent that any two black 
points inside hexagons surrounding these sites, are con- 
nected via black color. Thus, it terms of connectivity over 
black, two neighboring sites with Vq < s are completely 
similar to two neighboring black hexagons. Similarly, as 
can be seen from Fig. [2l for two neighboring sites with 
Vq < s and Vq > s, the centers of surrounding hexagons 
are disconnected. The same is true for two neighboring 
hexagons of opposite colors in the percolation problem. 
To complete the mapping, we note that, in percolation 
problem, the connectivity of two neighboring hexagons 
depends entirely on their colors, i.e., it does not depend 
on the color of the other neighbors. In the same way, in 
the problem of equipotentials, whether or not two neigh- 
boring sites are connected is determined exclusively by 
the signs oi Vq — e in these sites. 

APPENDIX B 

The form of the 5-matrix Eq. (I13p can be established 
with the help of Fig. [H Matrix S relates the incident, 
«i, i2, «c, and outgoing, oi, 02, Oc, amplitudes via trans- 
mission coefficients, xi, X2, and Xc. The form Eq. 
follows from the system of six equations, which include 
also the amplitudes oi, 02, and ac between the points 
of close contact of corresponding equipotentials. As seen 
from Fig. [Hthe amplitudes ai and a2 are related via xi 
as 

ai = xia2e'-'l'^ + yiii, 

02 ^ -~yia2e''l'^ + xiii. (Bl) 

The amplitudes Oc and ai are related via Xc 

Qc = Xc aie"^i + ycic, 
oi^-y^aie''!'' +Xcic- (B2) 
Finally, the amplitudes 02 and Uc are related via X2 

a2 = X2ace"^'' + ?;2«2, 

Oc = -y2ac e''''" + X2i2- (B3) 

In Eqs. (|B1[) - (|B3[) . the phases (f>i, (t>2, and (pc are the 
Aharonov-Bohm phases accumulated, respectively, by 
waves ai, 02, ac between the points of closed contact. 
Excluding ai, a2, ac from Eqs. (|Bip - (|B3p we recover 5*- 
matrix Eq. (I13p . in which X3 stands for Xc and 03 for 

4>c- 
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into the form 



FIG. 8: (Color online) A slice of a network Fig. |3]of width 
M = 8 is shown. Three amplitudes on the links to left of 
green dot are related to three amplitudes to the right of green 
dot via matrix X. The upper and the lower boundaries of the 
slice are connected by dashed-dotted lines manifesting that 
the amplitudes on these boundaries are the same by virtue of 
cyclic boundary conditions. Upon passing the white stripe, 
the vector of the amplitudes {Zi} is multiplied by the matrix 
Li. Upon passing the blue stripe the vector {Zi} is multiplied 
by L2. 



In order to derive Eq. we set in Eqs. (|BT |1 - (|B31) 
i2 — 02e"^^, as enforced by a contraction in Fig. [5^. 
Upon setting ii = 0, we find the proportionality coef- 
fcient between oi and ic- This recovers Eq. (jl6p . in 
which f = ip2 + <p2- Similarly, Eq. (|17p is recovered upon 
setting ii — 026^^, as shown in Fig. and relating the 
amplitudes oi and ic- 



APPENDIX C 

In this Appendix we illustrate the transfer-matrix 
method in application to the network model on trian- 
gular lattice. Consider, by analogy to Ref. H, a slice of 
length, N, shown in Fig. [S] It contains M — 8 inci- 
dent links, zl^\ and 8 outgoing finks, z'>^\ In general, 
the number of links must be M = 4m, where m is inte- 
ger. The slice shown in Fig. [5] can be obtained from the 
general network Fig. [3] by two vertical cuts through the 



(0) 



and Z^^"^ in 



centers of triangles. The bottom links Z^ 
Fig. [5] are connected to the corresponding top finks by 
dashed lines, reflecting the fact that these links must be 
identified with each other in order to impose a periodic 
boundary condition^. Scattering of waves at each node is 
described by the matrix S = S^,, Eq. ([7]), which relates 
the amplitudes (oi, 02,03) to (ii, 12, is). To adapt S — Si^ 
to the transfer-matrix algorithm, one has to recast Eq.fT^ 



02 ^ 










Ol 


03 J 







(CI) 



which connects the amplitudes to the left and to the right 
from the node. Direct calculation yields the following 

of S^. 



form of matrix X in terms of elements, stj 
( 



X 



V 



S21 



S31 



S12 



S12 



S12 



S22 
S12 



1 

S12 



S32 
S12 



S23 



S33 



S12 



Sl3 
S12 



S32S13 
S12 



(C2) 



For energies close to the critical e = 0, the 5-matrix is 
given by Eq. Q; S = St,. Making use of Eq. dZl), it is 
straightforward to find the energy dependence of X 



/2(l+e) -2{l + e) 1 \ 

2(1 +e) -3 2(1 -e) 
V 1 -2(1 -e) 2(1 -e)/ 



(C3) 



Operators Li and L2 act in "white" and "blue" stripes, 
respectively. Operator Li performs the transformation 
of the vector of amplitudes, {z\^^}, into {z[^^^^}, while 

L2 performs the transformation of {z["'^} into {Z^"-*}. 
Matrix forms of Li and L2 in terms of elements of matrix, 
X, are the following 



^ Xii X12 Xi3 

a;2i a;22 X23 

2^31 X32 X33 

1 



Li 



Xii X12 Xi3 

a;2i X22 X23 
2:31 a;32 2:33 



(C4) 





V • 














1/ 




^ X33 













3^31 









1 






















Xll 


a;i2 


Xl3 









L2 = 






X21 


X22 


X23 















X31 


X32 


X33 
























1 








2^13 













Xll 


a;i2 




\ a;23 













X21 


2^22 / 



(C5) 



with dots standing for zeroes. Specific form of L2 ac- 
counts for the cyclic boundary conditions in the vertical 
direction. 

In addition to the scattering at the nodes (green dots in 
Fig. [H]) , propagation along the network finks accumulates 
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random phases. If we denote the phases on the Unks 
crossing vertical hnes 2n by and those crossing hnes 

2n + 1 by V'i"'' : then all the random phases at the links 
can be taken into account by introducing diagonal phase 
matrices 

P^"^ = diagje*'^*"' , • • • , e'-^s"' } (C6) 

and 

P^") = diag{e*^i"' , • • • , e^'^'s"' }. (C7) 



Finally, the transfer matrix, T, of the slice Fig.[8]is given 
by the product 

T= f[ L2P("^LiPl"\ (C8) 

n=N-l 

which runs in the reverse order. This product is com- 
pletely analogous to the transfer matrix of the slice in 
the Chalker-Coddington model. 
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